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Abstract It is believed that the observed diffuse gamma 
ray emission from the galactic plane is the result of in- 
teractions between cosmic rays and the interstellar gas. 
Such emission can be amplified if cosmic rays penetrate 
, into dense molecular clouds. The propagation of cosmic 
' rays inside a molecular cloud has been studied assum- 
ing an arbitrary energy and space dependent diffusion 
coefficient. If the diffusion coefficient inside the cloud is 
significantly smaller compared to the average one derived 
for the galactic disk, the observed gamma ray spectrum 
appears harder than the cosmic ray spectrum, mainly 
due to the slower penetration of the low energy particles 
towards the core of the cloud. This may produce a great 
variety of gamma ray spectra. 
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1 Introduction 

The observed diffuse gamma ray emission from the galac- 
tic plane is believed to be the result of the decay of neu- 
tral pions produced during inelastic collisions of cosmic 
rays with interstellar gas. If, as a first order approxima- 
tion, one assumes that the cosmic ray spectrum is the 
same everywhere in the Galaxy, then the gamma ray 
emission is expected to be simply proportional to the 
gas column density. As a consequence, if the line of sight 
is intersecting regions of enhanced density, such as gi- 
ant molecular clouds, also the gamma ray emission is 
expected to be correspondingly enhanced. 

The importance of the detection of molecular clouds 
in gamma rays is widely recognize, especially in relation 
to the problem of the origin of cosmic rays. Molecular 
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clouds located in the vicinity of cosmic ray accelerators 
could provide a dense target for cosmic rays interactions, 
amplifying the resulting gamma ray emission and making 
easier the identification of cosmic ray sources mSEE]- 

On the other hand, even in the absence of an acceler- 
ator, molecular clouds embedded in the "sea" of galactic 
cosmic rays are expected to emit gamma rays. If cos- 
mic rays can freely penetrate the clouds, the high energy 
gamma ray spectrum is expected to mimic the slope of 
the cosmic ray spectrum and the total gamma ray lumi- 
nosity depends only on the cloud total mass. For this rea- 
son, molecular clouds can be used to probe the cosmic ray 
energy spectrum and its absolute flux in different parts 
of the Galaxy 5 3 4 . The condition for the detectability 
of these passive clouds with EGRET is M 5 /dk pc > 10, 
where M5 is the cloud mass in units of 10 5 M Q and dk pc is 
the source distance in kpc 0]. Since there are only a few 
clouds in the galaxy with such a large value, a more sen- 
sitive instrument like GLAST is needed for these studies. 

Of course, what discussed above is valid only un- 
der the assumption that cosmic rays freely penetrate 
the clouds. The issue of the penetration or exclusion of 
cosmic rays from clouds has been investigated in sev- 
eral papers [HltZIIHIj m which quite different conclusions 
have been drawn, going from the almost-free penetra- 
tion to the exclusion of cosmic rays up to tens of GeV. 
Since a theoretical determination of the cosmic ray dif- 
fusion coefficient is a very difficult task, we adopt here 
a more phenomenological approach: the diffusion coef- 
ficient is parametrized, and observable quantities (such 
as the gamma ray flux) capable of constraining it are 
proposed. 

In this paper we consider only passive clouds embed- 
ded in the diffuse galactic cosmic ray flux and located far 
away from cosmic ray accelerator. The case of clouds in- 
teracting with accelerators will be considered in a forth- 
coming paper. 
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2 Charcateristic time scales of the problem 

Before solving the equation that describes the transport 
of cosmic rays inside a magnetized, dense cloud, it is 
worth giving an estimate of the typical time scales in- 
volved in the problem. Consider a giant molecular cloud 
of radius R c i ~ 20 pc, mass M c ; ~ 2 10 5 M Q and average 
magnetic field B c i ~ 10 fiG. The average density (Hy- 
drogen atoms) of such cloud is thus n gas ~ 300 cm -3 , 
which gives a dynamical (free-fall) time for the system 
of the order of: 



Tdyn ~ {Gp)~ 



5.510' 



(300^) yr (1) 



where G is the gravitational constant and p the mass 
density. This free-fall time has probably to be considered 
as a strict lower limit to the lifetime of the cloud, since 
additional pressure support from fluid turbulence and 
magnetic field may inhibit the collapse |10j. 

The effectiveness of the cosmic ray penetration into 
the cloud depends on the interplay of several physical 
processes: (i) diffusion in the cloud magnetic field, (ii) 
advection due to turbulent bulk motion inside the cloud, 
(iii) energy losses in the dense cloud medium. Moreover, 
the cosmic ray density can be enhanced if a cosmic ray 
accelerator is embedded in the cloud [3], or if cosmic 
rays coming from outside the cloud are reaccelerated via 
Fermi-like processes that may take place in the magne- 
tized cloud turbulence 12101 ■ In the following we consider 
the cloud as a passive target for galactic cosmic rays and 
neglect any effect related to the possible presence of ac- 
celeration and/or reacceleration of cosmic rays inside the 
cloud. 

We parametrized the diffusion coefficient for protons 
of energy E propagating in the cloud magnetic field B 
in the following way: 



(2) 



where Dq = 3 10 27 cm 2 /s and <5 = 0.5 are the typical 
galactic values |TJ and the parameter x < 1 accounts for 
a possible suppression of the diffusion coefficient inside 
the turbulent cloud medium. In general, the values of \ 
and S will depend on the power spectrum of the magnetic 
field turbulence. For such a choice of parameters, one can 
estimate the proton diffusion time, namely, the time it 
takes a proton to penetrate into the core of the cloud: 
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To study the effective propagation of cosmic rays into 
cluods, it is instructive to compare the diffusion time 
with the energy loss time. In the dense cloud environ- 
ment, cosmic ray protons suffer energy losses due to 



ionization and nuclear p-p interactions. Above the en- 
ergy threshold for pion production Eth ~ 300MeV nu- 
clear interactions dominate. Since both the cross section 
tjpp ~ 40 mb and inelasticity k ~ 0.45 of this process are 
not changing significantly over a broad range of proton 
energies from ~ 1 GeV to hundreds of TeVs, the proton 
lifetime is almost energy independent: 
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Cosmic rays can also be transported by the fluid 
turbulence which is known to be present in molecular 
clouds. Several molecular lines are observed in the direc- 
tion of clouds, and their width Av reflects the velocity 
of internal turbulent motions. The line width is known 
to correlate with the cloud size according to the relation: 
Av oc i?^ 5 ^U], which for the cloud sizes considered here 
provides a velocity of a few km/s. The time scale for this 
advective transport can be roughly estimated as: 



Tadv ^ 4 1U ^ 2()pc y ^ 5km/s 



/ Av 



yr 



(5) 



The energy dependence of all the time scales con- 
sidered above is shown in Figure ^ (left panel) , where 
the solid and dashed thick lines represent the diffusion 
time (with \ = 1 and \ — 0-01 respectively), while the 
horizontal lines refers to the dynamical, advective and 
energy loss time scales (from top to bottom). Several 
comments are in order. First of all, the dynamical life- 
time of the cloud is the longest time scale for all the rel- 
evant energies. This means that it is possible to search 
for a steady state solution of the problem. Second, the 
advection time is comparable with the dynamical time, 
but it is always significantly longer than both the diffu- 
sion and the energy loss time scales. Thus, we can safely 
ignore the advection term in the cosmic rays transport 
equation. As a consequence, the degree of penetration 
of cosmic rays inside the cloud can be roughly estimated 
by comparing the diffusion and energy loss times. To this 
purpose, in Figure ^ (left panel) two different values of 
the parameter \ are considered. If x = 1 the diffusion 
time is shorter than the energy loss time at all the con- 
sidered energies. Thus, if the diffusion coefficient inside 
the cloud is not suppressed with respect to the galactic 
value, cosmic rays can easily penetrate into the cloud. 
On the other hand, if diffusion inside the cloud is sig- 
nificantly suppressed (x = 0.01), the energy loss time 
becomes shorter than the diffusion time at energies be- 
low E* ~ 10 — lOOGeV. This means that only cosmic 
rays with energy above E* can penetrate into the cloud 
before losing their energy. This is a very improtant fact, 
since as we will demonstrate in the following sections, 
the exclusion of low energy cosmic rays plays a crucial 
role in shaping the gamma ray spectrum of molecular 
clouds. 

Finally, it is worth stressing that the value of E* in- 
creases if one considers a realistic density profile for the 
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Fig. 1 LEFT PANEL: typical time scales for cosmic rays in a giant molecular cloud with average number density ~ 
300 cm -3 , radius 20 pc and average magnetic field 10 ^iG. Solid and dashed lines represent the diffusion time, horizontal 
lines represent the dinamical, advective and energy loss times (top to bottom). RIGHT PANEL: cosmic ray spectrum in 
the cloud center for different values of the parameter \, describing the suppression of the diffusion coefficient with respect 
to the measured galactic value (see text for details). 



cloud instead of average quantities. Despite the fact that 
density profiles cannot be easily extracted from available 
observations, it is well known that molecular clouds con- 
tain cores of size ~ 1 pc or less, in which density can 
reach very high values n gas ~ 10 5 cm~ 3 . In this case, 
a proton which is approaching the center of the cloud 
meets a denser and denser environment. The value of the 
magnetic field is also increasing towards the cloud cen- 
ter, reaching in the densest regions values of hundreds of 
\xG or even more These facts make the energy loss 
time shorter and the diffusion time longer, leading to a 
more efficient exclusion of cosmic rays from clouds cores. 



3 Solution of the transport equation 

As shown above, the equation describing the cosmic ray 
transport in a cloud is the steady state diffusion-losses 
equation, which in spherical simmetry reads: 



r^or 



D(R, E)R 



2 d£ 
dR 



d_ 
~dE 



E(R)f 







(G) 



where f(R, E) is the space and energy dependent particle 
distribution function, D(R,E) is the diffusion coefficient 
parametrized as described in the previous section and 
E(R,E) = dE/dt represents energy losses. 

We parametrize the density profile as follows: 



n H (R) 



(7) 



where n is the central density and R c the core radius, 
assumed to be 1/2 parsec. This is of course a simpli- 
fied assumption, since molecular clouds can exhibit very 
irregular density profiles. However, we will show in the 
following that the mechanism of cosmic ray exclusion 
from clouds can work for very different density profiles, 
from flat (a — 0) to very peaked ones (a ^> 0). 

In calculating the energy loss term we consider only 
inelastic proton-proton collisions, since this is the dom- 
inant process above the threshold for the production of 
pions (E th ~ 280 MeV). The loss term in Eq. de- 
pends on the density profile through Eq. Q , and thus it 
is space dependent. 

For the magnetic field profile we use the results from 
|12|. in which Zeeman measurements of magnetic field 
strength in molecular cloud cores are reported. A corre- 
lation between magnetic field strength and gas density 
is observed and can roughly be fitted by: 



B ~ 100 



( nH ) 
V10 4 cm" 3 / 



1/2 



(8) 



Though this correlation has been found for molecular 
cloud cores with density exceeding nu ~ 10 3 cm -3 , it 
provides reasonable values also for low density regions 
(tens of \xG for typical average cloud densities of a few 
100cto~ 3 ). Thus, we assume that Eq. (JSJ is valid in the 
entire density interval. It is worth noticing that in our 
model the diffusion coefficient is space-dependent, since 
it depends on the magnetic field as given by Eq. (J2J. 

We solved Eq. © numerically, using an implicit scheme 
and assuming, as boundary condition, that the cosmic 
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ray spectrum outside the cloud must match the galactic 
cosmic ray spectrum |H] . The galactic spectrum is in turn 
assumed to be equal to the locally observed cosmic ray 
flux: 

Jg*(E) = 2.2 ^— J cm^s-hr-'GeV- 1 (9) 

The effective exclusion of cosmic rays from the cloud 
cores is demonstrated in Fig. ^ (right panel) , where the 
cosmic ray spectrum in the cloud center is plotted for 
different values of the parameter \ (here S is set equal to 
the galactic value 0.5). In obtaining the result, a cloud 
with mass 2 x 10 5 M Q , radius 20 pc and a flat density 
profile (this correspond to a spatially constant magnetic 
field of ~ 15/iG) has been considered. If the diffusion co- 
efficient inside the cloud is not suppressed with respect to 
the galactic value (x — 1), then the cosmic ray spectrum 
in the cloud center is basically indistinguishable from the 
galactic cosmic ray spectrum. On the other hand, is dif- 
fusion is significantly suppressed (x << 1), cosmic rays 
with energy above - 10 GeV (x = 0.01) or ~ 100 GeV 
(x = 0.001) cannot penetrate the cloud. As we will show 
in the following, this fact has important implications for 
the estimate of the spectrum and intensity of the gamma 
ray emission expected from molecular clouds. 

Cosmic ray protons propagating inside a molecular 
cloud also produce secondary electrons during inelastic 
interactions in the intercloud medium. These electrons 
contribute to the overall gamma ray emission of the cloud 
via Bremsstrahlung. Once the steady state proton spec- 
trum has been obtained, we calculated the injection spec- 
trum Q e (E e ) of the secondary electrons by using the 
analytical fits provided in ^21- The steady state spec- 
trum of secondary electrons can be obtained using again 
Eq. © appropriately modified as follows: i) the injec- 
tion term Q e (E e ) must be added on the left side; ii) 
the loss term E e (R, E) is now dominated by Coulomb, 
Bremsstrahlung and synchrotron losses and it is both 
space and energy dependent; Hi) as boundary condition 
we assumed that the particle distribution function for 
secondary electron vanishes outside the cloud, where the 
gas density is low with inefficient production of secon- 
daries. 



4 Gamma ray spectra 

Fig. [21 (left panel) shows the gamma ray spectra for a 
cloud of mass M — 2 x 1O 5 M and radius R c i — 20pc. 
A flat density profile is assumed for a cloud at a dis- 
tance of 1 kpc. The thick lines represent the contribu- 
tion to the total gamma ray emission from 7r°-decay 
(calculated following JS]), while the thin lines represents 
the Bremsstrahlung contribution. Calculations have been 
performed adopting a diffusion coefficient with S = 0.5, 
but assuming three different values for \ = 1, 0.01, 0.001 
(curves top to bottom respectively). The Bremsstrahlung 



contributions become significant only below photon en- 
ergies of - 100.MW. 

The solid line represents the situation in which cos- 
mic rays can freely penetrate the cloud, namely x = 1 
(no suppression of the diffusion coefficient). This gives 
the maximum possible gamma ray luminosity for a pas- 
sive cloud immersed in the galactic cosmic ray sea. If 
the diffusion coefficient is suppressed (dashed and dotted 
lines, corresponding to \ = 0.01, 0.001) the total gamma 
ray luminosity is reduced, especially at low energies (~ 1 
GeV), while at high energies the canonical spectrum for 
a passive cloud filled by cosmic rays is recovered. This 
reflects the fact that high energy cosmic rays can freely 
penetrate inside clouds. At ~ 1 GeV, the suppression 
of the gamma ray flux is roughly a factor of ~ 2 for 
X = 0.01 and a factor of ~ 5 for x — 0.001. Moreover, 
also the shape of the spectrum in modified, appearing 
flatter for lower values of \ (stronger suppression of dif- 
fusion) . This fact may have very important implications 
for future GLAST observations of molecular clouds. In 
particular, the observation of gamma ray spectra harder 
than the ones expected from an isolated (passive) cloud 
pervaded by galactic cosmic rays can be interpreted in 
two different ways: i) the galactic cosmic ray spectrum 
at the location of the cloud is different than the one 
measured locally; ii) propagation effects inhibit the pen- 
etration of low energy cosmic rays in the cloud, making 
the resulting gamma ray spectrum harder. 

Molecular clouds constitute also potential targets for 
Cherenkov telescopes arrays, operating at photon ener- 
gies greater than 100 GeV. The apparent angular size of 
a cloud with radius R c i ~ 20 pc located at d = 1 kpc is 
■d c i ~ 2R c i/d ~ 2.2°. This is significantly smaller than 
the telescope field of view (e. g. ~ 5° for HESS) and 
much larger than its angular resolution (~ 0.1°), thus 
Cherenkov telescopes can effectively map the gamma ray 
emission from clouds. Notably, a core with radius of half 
parsec will subtend an angle comparable with the tele- 
scope angular resolution. For this choice of the param- 
eters, the gamma ray flux (only 7r°-decay contribution) 
for a cloud of mass 2 x 1O 5 M is shown in Fig. [21 (right 
panel). The diffusion coefficient is the galactic one sup- 
pressed by a factor of 100 and the parameter a describ- 
ing the slope of the density profile is varied. In the top 
panel the emission from the whole cloud is shown. Thick 
lines represent the spectrum for different values of a: 
(solid), 1 (dashed) and 2 (dotted). The thin line rep- 
resents the spectrum one would observe if cosmic rays 
could freely penetrate the cloud. It can be seen that the 
emission from the whole cloud is not depending strongly 
on the density profile. The situation is much different if 
one considers the radiation from the inner parsec region 
(Fig. [21 bottom panel), which in this particular case cor- 
respond with the radiation received within one angular 
resolution of the Cherenkov telescope. In this case the 
spectrum is strongly dependent on the assumption made 
on the density profile, especially at very high energies 
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Fig. 2 LEFT PANEL: gamma ray emission from a cloud with M — 2 x 10 5 Mq, R ct — 20 pc and a flat density profile. 
The cloud distance is 1 kpc. Thick lines: 7r°-decay gamma rays, thin lines: Bremsstrahlung gamma rays. RIGHT PANEL: 
Thick lines: gamma ray emission (only 7r°-decay component) from clouds with different density profiles and \ = 0.01. Thin 
lines show the gamma ray emission one would observe if cosmic rays could freely penetrate the cloud. 



above ~ 10 GeV. This is because, for peaked density 
profiles [a > 0) the exclusion of cosmic rays from cloud 
cores is much more effective with respect to the case of 
a flat profile (a = 0), due to the enhanced energy losses 
and reduced diffusion there. This makes the gamma ray 
spectrum harder. On the other hand, since the cores are 
very dense, the suppression of the cosmic ray density 
is compensated by the higher efficiency of production 
of gamma rays. This explains why a higher level of the 
gamma ray emission from the densest cores is predicted. 

Thus, we arrive at the important conclusion that the 
shape of the gamma ray spectrum from molecular cloud 
cores is determined not only by the diffusion properties 
of cosmic rays but also by the actual shape of the cloud 
density profile. 

This statement is confirmed in Fig. [3 where the slope 
of the gamma ray spectrum at photon energies of 100 
GeV (top panels) and 1 TeV (bottom panels) is plot- 
ted for different combinations of model parameters. In 
the left panel we fix 6 = 0.5 but we change the slope 
of the density profile a, while in the right panel, a is 
kept fixed and the energy dependence of the diffusion 
cocfficiocn 6 is varied. In both plots x — 0.01. Triangles 
refer to the spectrum from the whole cloud and circles 
refer to the core. Again, we emphasize that, although the 
total gamma ray spectrum from the whole cloud is not 
strongly affected by the choice of the model parameters, 
the cloud cores may show a great variety of energy spec- 
tra. As a consequence, gamma ray observations in the 
TeV energy domain can be very useful to constrain pa- 
rameters. In principle, if the cloud density profile could 
be extracted from independent observations (e.g. CO line 



emission) the shape of the gamma ray spectrum would 
tell us important information about the diffusion coeffi- 
cient. Of course, the determination of the density profile 
is a extremely difficult task, since clouds show very com- 
plex structures. 



5 Conclusions 

In this paper we considered a giant molecular cloud em- 
bedded in the diffuse galactic cosmic ray flux. Assum- 
ing that the cosmic ray propagation inside the cloud 
proceeds in the diffusive regime, we studied the exclu- 
sion/penetration of cosmic rays into molecular clouds. 
Results can be summarized as follows: 

• if the diffusion coefficient inside the cloud is equal to 
the measured galactic one, cosmic rays can freely pen- 
etrate the cloud and the resulting high energy gamma 
ray emission has a spectral shape that closely resem- 
bles the one of the galactic cosmic rays. 

• the gamma ray emission above 100 MeV is dom- 
inated bi 7r°-decay. Bremsstrahlung emission from 
secondary electrons is relevant only at smaller en- 
ergies. 

• if the diffusion coefficient is suppressed with respect 
to the galactic one, cosmic rays can be effectively ex- 
cluded from clouds. In particular, for a suppression of 
the diffusion coefficient of a factor of ~ 100, the ex- 
clusion becomes relevant at energies of tens-hundreds 
of GeV. 
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Fig. 3 Slope of the gamma ray spectrum for different sets of parameters. Triangles refer to the spectrum of the whole cloud, 
circles refer to the spectrum observed by looking at the cloud core. 



• the exclusion of cosmic rays from clouds results in a 
suppression of the gamma ray flux, especially at ~ 
GeV energies. This can have important consequences 
for the forthcoming GLAST observations. 

• Cherenkov telescopes such as HESS or VERITAS has 
the capability to map the gamma ray emission from 
clouds and to resolve the pc-scale cores. In particular, 
the shape of the gamma ray spectra from cloud cores 
may strongly depend on both the diffusion properties 
of cosmic rays and the shape of the cloud density 
profile. The effect of the density profile in shaping the 
gamma ray spectrum of molecular clouds was never 
considered before. 

This work has been stimulated by the recent detec- 
tion of some molecular clouds in the galactic center re- 
gion by HESS [El and by the forthcoming observations 
of the galactic plane by GLAST. It is beyond any doubt 
that future observations of molecular clouds by these two 
instruments, covering a broad energy interval from GeV 
to multi TeV photon energies, will provide deep insight 
into the problem of the origin of galactic cosmic rays. 
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